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Abstract 

We investigate the nonextensive g-distribution function for a gas in the pres- 
ence of an external field of force possessing a potential U(r). In the case of 
a dilute gas, we show that the power law distribution including the potential 
energy factor term can rigorously be deduced based on kinetic theoretical ar- 
guments. This result is significant as a preliminary to the discussion of long 
range interactions according to nonextensive thermostatistics and the underly- 
ing kinetic theory. As an application, the historical problem of the unbounded 
isothermal planetary atmospheres is rediscussed. It is found that the maxi- 
mum height for the equilibrium atmosphere is z m ax = KbT /mg(l — q). In the 
extensive limit, the exponential Boltzmann factor is recovered and the length 
of the atmosphere becomes infinite. 
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In gas dynamics it is very important to investigate how the molecular motion is modified 
by force-fields different from those exerted by the containing vessel or even by the other 
particles of the gas. Particular examples are ions in an external magnetic field and a gas in 
the earth's gravitational field [1,2] 

As widely known, a classical gas under steady state conditions and immersed in a con- 
servative force field, F = — VC/(r), is described by a distribution function that differs from 
the Maxwellian velocity distribution by an extra exponential factor involving the potential 
energy In this case, the total equilibrium distribution function reads 

where m is the mass of the particles, T is the temperature and n is the particle number 
density in the absence of the external force field. In addition, since this distribution function 
is normalized, it is easy to see that the number density is given by 



n(r) = no exp 



(2) 



k B T_ 

where the factor, exp[— (/(rj/^T], which is responsible for the inhomogeneity of f(r,v), is 

usually called the Boltzmann factor. Expression (1) follows naturally from an integration of 

the collisionless Boltzmann's equation 

df df F df 
at Or m <7v 

when stationary conditions are adopted along with the assumption that the total distribution 
can be factored 

/(M) = /o(#), (4) 

where fo(v) represents the Maxwell equilibrium distribution function, and x( r ) is a scalar 
function of r. As one may show, after a simple normalization, the resulting expression for 
x( r ) is exactly the Boltzmann factor for the potential energy of the external field, namely: 



X(r) = exp 



(5) 



2 



and combining this result with equation (4) the Boltzmann stationary distribution (1) is 
readily obtained. 

On the other hand, recent efforts on the kinetic foundations of the g-nonextensive statis- 
tics proposed by Tsallis [3,4] lead to an equilibrium velocity distribution of the form [5,6] 



In this expression the g-parameter quantifies the nonadditivity property of the associated gas 
entropy whose main effect at the level of the distribution function is to replace the standard 
Gaussian form by a power law. The quantity B q is a g-dependent normalization constant 
whose expression for 1/3 < q < 1 is given by 



which reduces to the Maxwellian result in the limit q — 1. As explained in Ref. [5], the above 
power law distribution can be deduced from two simple requirements: (i) isotropy of the 
velocity space, and (ii) a suitable nonextensive generalization of the Maxwell factorizability 
condition [7], or equivalently, the assumption that F(v) = f(v x )f(v y )f(v z ). For q > 0, the 
above distribution function satisfies a generalized H-theorem, and its reverse has also been 
proved, that is, the collisional equilibrium is given by the Tsallis' q-nonextensive velocity 
distribution [6]. 

In the last few years, several applications of this equilibrium power law velocity distri- 
bution have been done in many disparate branches of physics [8-14]. It is worth notice that 
the widely spread belief that Tsallis equilibrium distribution cannot be applied to Hamil- 
tonian systems seems to be a profound misunderstanding on the foundations of statistical 
mechanics and kinetic theory. Indeed, the BG canonical emsemble approach is valid only 
for sufficiently short-range interactions. It fails when gravitational or unscreened Coulom- 
bian fields are present, that is, to the class of systems where the usually assumed entropy 
additivity postulate is not valid. In other words, when the dynamics plays a nontrivial role, 
the system does not relaxes to the BG distribution, but they evolve to a non-Gaussian ve- 



mv 
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fo(v) = B q l-(l-q) 



2k B T 



(6) 
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locity distribution which can be fitted by the monoparametric class of Tsallis' nonextensive 
statistics (for a more detailed discussion see [15,16]). 

In this letter we investigate how the potential energy term can rigorously be introduced 
in the g-distribution with basis on the kinetic approach. More precisely, it is shown that 
an analytical expression for the equilibrium distribution of a dilute gas under the action 
of a conservative field may also be calculated as a stationary solution of the collisionless 
Boltzmann equation. This result is significant as a preliminary to the discussion of long 
range interactions according to nonextensive thermostatistics and the underlying kinetic 
theory. As an illustration we discuss the case of planetary atmospheres. It is found that the 
problem of an infinite isothermal atmosphere is naturally solved in this extended framework. 
For a given value of q < 1, the extension of the atmosphere becomes finite and its length is 
uniquely determined by the temperature and mass of the molecular components. 

Let us now consider a spatially inhomogeneous dilute gas supposed in equilibrium at tem- 
perature T. It is immersed in a conservative external field in such a way that f(r,v)d 3 vd 3 r 
is the number of particles with velocity lying within a volume element d 3 v about v and 
positions lying within a volume element d 3 r around r. In this case, we see from (7) that the 
stationary Boltzmann equation can be rewritten as 



In order to introduce the nonextensive effects we first recall that the factorizability 
condition cannot be applied in this extended framework. This means that the standard 
starting assumption (see equation (4)) must be somewhat modified. In the spirit of the 
g-nonextensive formalism, a consistent g-generalization of (4) is 



where B q is the constant expression given by (7) which has been introduced for mathematical 
convenience, and the functions g-exp and g-log, e q (f), ln q (f), are defined by 




0. 



(8) 




(9) 




(10) 
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Note that in the limit q — > 1 the above identities reproduce the usual properties of the 
exponential and logarithm functions. In addition, e q (\n q f) = \n q (e q (f)) = f, and as should 
be expected from (9), the factored decomposition (4) is readily recovered in the extensive 
limit. In what follows, the properties of g-exponential and g-log differentiation 

= /-.*, (12) 

will also be extensively used. In particular, for the total g-distribution (9), we obtain 

TYUV 

V„/ = -e 3 W^. (14) 

Now, substituting V r / and W v f into the stationary Boltzmann equation (8), and performing 
the elementary calculations one obtains 

X - q V X .dr = --1-VU.dr (15) 
k B T 

the solution of which is 

* w ^«(-Sf + c )- (16) 

where C is an arbitrary constant. 

Now, inserting (16) into (9), and integrating the result in the velocity space it follows 
that 



/ 



B n e n 



d s v = n(r). (17) 



Finally, by substituting the expression of f and considering a region where U(r) = 0, one 
finds 



B q je q [-^ + C\d A v = n 0) (18) 



r 2 

2k bT 



and from normalization condition, no — J fo(v)d 3 v, it follows that the unique allowed value 
for the integration constant is C = 0. Consequently, (16) becomes 

u( r y 



X(r) 



k B T 



(19) 



which is the g-generalized Boltzmann factor. 

Finally, inserting this result into (9), we obtain the complete g-distribution function in 
the presence of an external field 

V(l-9) 



f(r,v) = B q 



!-(!-?) 



mv 



U(v) 



= B q e q {-E/k B T), 



(20) 



2k B T k B T 

where E is the total energy of the particles. It thus follows that a generalized g-exp factor 
for Tsallis' nonextensive thermostatistics can exactly be deduced if the standard approach is 
slightly modified. The particle number density, n(r) = / /(r, v)d 3 v , also becomes a function 
of the position given by 

(5-3g)/2(l-g) 



n(r) = n 



1 -(!-?) 



knT 



= n [e q (-U/k B T)]^ 



(21) 



and as should be expected, in the extensive limit (q — > 1), the standard exponential 
Maxwellian expressions for /(r, v) and n(r) are readily recovered. 

At this point it is interesting to analyze some application of the above result. We recall 
that the most popular problem of a gas in a force-field is the planetary atmosphere. In the 
standard simplified treatment, the temperature is uniform and the tree-dimensional motion 
occurs under the action of a constant gravitational field g along the z-direction. Thus, each 
molecule has a potential energy, U(z) = mgz, so that the concentration of particles, or 
equivalently, the gas density (p = nm) is given by the familiar barometric formula (see Eq. 



(2))- 



p(z) = po exp 



mgz 



k B T 



(22) 



with the pressure obeying a similar exponential law in virtue of the standard equation of 
state for a perfect gas. The simplest conclusion based on the above expression is that an 
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ideal isothermal atmosphere would extend to indefinite distances from the planet surface. 
This result disagree completely with the existing observations even by considering that the 
upper-most part of the earth atmosphere is nearly in thermodynamic equilibrium. Since 
there is no a natural boundary in this Boltzmann like picture, a possible solution to this 
problem is to assume that the atmospheres of the planets are actually finite, although being 
only in a pseudo-equilibrium state with the particles continuously streaming off into space 
at a very slow rate (evaporation). 

How does this basic prediction is modified when an equilibrium g-distribution is assumed? 
The nonextensive answer to this question is immediate. The new distribution itself is given 
by (see (20)) 



while from (21), the matter density assumes the form 



1 1/(1-9) 



(23) 



P(z) = Po 



1 -(!-?) 



mgz 



(5-3g)/2(l-g) 



(24) 



For comparison, in Fig. 1 we plot the gas density as a function of z for some values of the q- 
parameter including the Boltzmannian case. Note that the standard exponential curve (q = 
1) is replaced by the characteristic power law behavior of Tsallis' nonextensive framework. 
In a more realistic treatment, as in the case of mixed atmospheres, the above distribution 

(23) -(25) may separately be applied to each component. As in the Boltzmann case, heavier 
gases are more concentrated near the surface of the earth than the lighter ones. In particular, 
this means that at greater heights the fraction of oxygen is smaller than helium as compared 
to the corresponding fractions at the surface. Such a result remains true as long as the 
maximum height allowed for a given component is not attained. In fact, for values of q 
smaller than unity, the positiviness of the argument appearing in the power law implies that 

(24) exhibits a "thermal cut-off" to the allowed values of the z-coordinate which defines the 
length of the atmosphere. It thus follows that the maximum height of the atmosphere (or 
more generally, the length of any isothermal layer) is 



We see that for given values of the temperature and q < 1, the maximum height associated 
to a specific component is inversely proportional to its molecular mass. In table 1, we display 
the predicted values of z max for a temperature T = 298, IhK and two different kind of gases 
(Oxygen and Hydrogen). In figure 1, one may see two different plots showing the finiteness 
of the atmosphere for values of q smaller than unity. Note that without limit 

when q — > 1 recovering the Maxwellian result. As one may check, the predicted atmosphere 
is infinite for all values of q > 1 since the matter density, as given by equation (24), does 
not vanish at any finite height. 

The above solution to the isothermal atmosphere with infinite length resembles a similar 
solution proposed in galactic dynamics to the collisionless stellar isothermal sphere. In 
the latter case, the Maxwell-Boltzmann distribution also predicts an infinite system whose 
integrated density profile yields an infinite total galactic mass. This theoretical problem 
was also solved with help of the nonextensive power law distribution [8]. In the case of 
planetary atmosphere, the nonextensive prediction of a finite length, as well as that heavier 
components are associated to a smaller maximum height could be tested by experiments 
using a mass spectrograph in atmospheric balloons. 



TABLES 

TABLE I. Limits to the height of the atmosphere 



Gas (T = 298, 15 K) q z max (Km) 

Hydrogen: 1.0 oo 

0.9 1239 

0.8 619.5 

0.7 413 

Oxygen: 1.0 oo 

0.9 77.4 
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FIGURES 




FIG. 1. Maximum height of the isothermal atmosphere. The gas density as a function of z for 
some selected values of the q parameter. In the left panel, the solid line (q = 1) is the exponential 
curve for an infinite atmosphere predicted by the standard Boltzmann distribution. For values of 
q < 1 all curves decrease rapidly toward the maximum value of z which is a physical cut-off present 
in the power law behavior. The right panel displays the logarithm of the gas density versus z 
making more clear the existence of this cut-off. Note that the standard curve (q = 1) is a straight 
line going to infinite. However, for values of q < 1, there is a knee causing an intersection of the 
curves at finite values of the z axis. 

As a conclusion, we call attention for some simple physical situations where the power law 
distribution (24) can play an interesting role. The first example is the question concerning 
the escape of molecules to infinity. This phenomenon seems to be more easily understood 
in this extended framework since the power law distribution naturally defines a limiting 
boundary frontier where the atmosphere terminates (see Table 1 and Fig. 1). We recall that 
in the standard treatment, the rate of loss to infinity is based on the assumption that the 
upper part of the atmosphere extends in isothermal equilibrium only until a certain height 
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h. In other words, unlike in the g-nonextensive framework, there is a cut-off introduced 
by hand in the infinity isothermal atmosphere predicted by the Boltzmann approach. A 
more quantitative analysis of evaporation taking into account the earth curvature will be 
presented in a forthcoming communication. 

Other examples of forces, like the ones involving giroscopic terms may also be trivially 
added to the power law distribution (24). For rotating frames with constant angular velocity, 
for instance, the effect is just to add the term —l/2muj 2 l 2 to the potential energy, where / is 
the perpendicular distance from the axis of rotation. As usual, this term simulates a slight 
change in the potential energy due to gravity [17]. 

Finally, we also remark that for colloidal particles in solutions, the gravity acts on the 
particles which are buoyed up by the liquid in which they are suspended. As a consequence, 
instead of the weight mg the net force is mg(p p — pi) / p p , where p p and pi are the particles 
and liquid densities, respectively. In this case, the resulting nonextensive expression for n(z) 
could experimentally be tested by repeating the experiments with colloidal suspensions as 
originally done by Perrin (see [18] and references there in). Naturally, the inclusion of other 
physical effects leading to independent probes of the q-statistics take even greater interest. 
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